#Aula 5

#Exercicios 5 - Criacao de graficos

#Exercicio 5.1 - Editando alguns parametros graficos

par(mfrow=c(2,2))    #Janela grafica com 2 linhas e 2 colunas, onde caber 4 graficos
par(las=1)           #Legenda do eixo na horizontal
par()$mar            #para saber quais sao as margens
par (mar=c(1,1,1,1))  #Arrumar o espaco entre os eixos
par(tcl=0.5)        #Mudar os ticks, colocaar para dentro
par()$cex           #Para saber qual o tamanho da fonte que estava
par(cex=1.5)        #Alteracao no tamanho de todas as fontes
par()$cex.axis
par(cex.axis=1.5)   #Alterao o tamanho da fonte dos eixos
par()$cex.lab
par(cex.lab=2)      #Alterar tamanho da fonte do titulo dos eixos
par()$cex.main    
par(cex.main=3)    #Alterar o tamanho do titulo
par()$lwd
par()$mai         # qual e o espaco entre o grafico e a area da figura
par(mai=c(3,3,3,3))
par(bty="l")       #formato da caixa do grafico. Como so quero o x e y, a caixa tem formado de L

# Lendo os dados
esa=read.csv("esaligna.csv", header = T, sep=",", dec=".", as.is=T)
str(esa)
head(esa)
summary(esa)

#Grafico
par(las=1, tcl=0.5,cex=0.8, mar=c(8,8,4,8),bty="l")
plot(esa$dap ~ esa$ht, xlab="DAP (cm)", ylab="Altura (cm)", main="Relao entre DAP e altura")

#Exercicio 5.2

# Lendo os dados
esa=read.csv("esaligna.csv", header = T, sep=",", dec=".", as.is=T)
str(esa)
head(esa)
summary(esa)

#Explicacao das funcoes utilizadas

locator()     #Usado para especificar um ponto onde quero colocar a legenda A e B para o grafico de boxplot e barplot, respectivamente
?mtext()
mtext("a", side=2, at=28.0)   #Inserir a letra a, no canto superior esquerdo (2), na posicao 28

#Janela grafica com 1 linhas e 2 colunas, para os 2 graficos
par(mfrow=c(1,2))

#Criacao do boxplot
par(las=1, tcl=0.5,cex=0.8,bty="l")
boxplot(esa$dap ~ esa$talhao, xlab="Talho", ylab="DAP (cm)")

#Colocando a letra A no grafico do boxplot
locator()
mtext("a", side=2, at=28.0)

##Dados para o segundo grafico:

#Calcular media do talhao
med.talhao = aggregate(esa$dap, list(esa$talhao), FUN=mean)
colnames(med.talhao) = (c("Talhao","Media DAP"))
med.talhao = as.matrix(med.talhao)
med.talhao
str(med.talhao)

#Calculando o desvio padrao do talhao
desvpad.talhao = aggregate(esa$dap, list(esa$talhao), FUN=sd)
colnames(desvpad.talhao) = c("Talhao", "desv.pad.DAP")
desvpad.talhao = as.matrix(desvpad.talhao)
str(desvpad.talhao)
desvpad.talhao

#Unindo as matrizes
mat.DAP = matrix(c(med.talhao, desvpad.talhao), ncol=4)
mat.DAP
colnames(mat.DAP) = c("Talhao", "Media", "Talhao", "Desv.Pad")
mat.DAP

# Criando o grafico 2 

par (mar=c(10,6,1,1))
plot(mat.DAP [,c(1:2,4)], xlab="Talho", ylab="DAP (cm)", ylim=range(0:25), pch=(19), yaxp="n")
axis(side=1, labels=c("16", "17", "18", "22", "23", "32"))

#Calculando os segmentos
desv.menor = (med.talhao) - (desvpad.talhao)
desv.menor
desv.maior = (med.talhao) + (desvpad.talhao)
desv.maior
segments(x=c(16,17,18,22,23,32), y=desv.menor[,2], x1=c(16,17,18,22,23,32), y1=desv.maior[,2])
















